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We calculate the unregularized monopole and dipole contributions to the self-force acting on a 
particle of small mass in a circular orbit around a Schwarzschild black hole. From a self-force point 
of view, these non-radiating modes are as important as the radiating modes with I > 2. In fact, we 
demonstrate how the dipole self-force contributes to the dynamics even at the Newtonian level. 

The self-acceleration of a particle is an inherently gauge-dependent concept, but the Lorenz gauge 
is often preferred because of its hyperbolic wave operator. Our results are in the Lorenz gauge and 
are also obtained in closed form, except for the even-parity dipole case where we formulate and 
implement a numerical approach. 



I. INTRODUCTION 

The capture of solar-mass compact objects by massive 
black holes residing in galactic centers has been identi- 
fied as one of the most promising sources of gravitational 
waves for the Laser Interferometer Space Antenna 
The need for accurate templates for signal detection and 
source identification is currently motivating an intense 
effort from many workers to determine the motion of a 
relativistic two-body system in the small mass-ratio limit, 
without relying on slow-motion or weak-field approxima- 
tions; for a review, see Ref. 0. The work presented in 
this paper is part of this larger effort. 



Gravitational self-force and the MiSaTaQuWa 
equations of motion 

Consider a small body of mass m in orbit around a 
much larger black hole of mass M. In the test-mass limit 
(to — > 0) the motion of the small body is known to fol- 
low a geodesic in the spacetime geometry of the large 
black hole 0, 0> 0, El • But as the mass of the smaller 
object is allowed to increase, deviations from geodesic 
motion become noticeable; these are associated with im- 
portant physical effects such as radiation reaction and 
finite-mass (conservative) corrections to the orbital mo- 
tion. In a sense, the motion is now geodesic in the per- 
turbed spacetime that contains both the black hole and 
the orbiting body. If g a g denotes the unperturbed metric 
of the central black hole, and if h a p denotes the pertur- 
bation produced by the orbiting body, then the motion 
is formally geodesic in the perturbed metric g a g + h a @. 
When viewed from the background spacetime, the motion 
appears to be accelerated, and the agent that produces 
the acceleration is a gravitational self-force acting on the 
particle. 

To turn these considerations into concrete equations of 
motion, it is desirable to formulate an approximation in 
which the details of the small body's internal structure 



have a negligible influence on the body's orbital motion. 
In this approximation the body is modeled as a point par- 
ticle possessing mass but no higher multipole moments, 
and its motion is fully described in terms of a world line 
7. But to formulate equations of motion for this world 
line becomes problematic, as the field h a g produced by a 
point particle necessarily diverges at the position of the 
particle. This means that an afiine connection cannot be 
defined on the world line, and that the statement "the 
particle follows a geodesic of the perturbed spacetime" 
does not make immediate sense. 

The task of regularizing h a g near the world line and 
formulating meaningful equations of motion for the point 
particle was undertaken by Mino, Sasaki, and Tanaka [jj, 
and also by Quinn and Wald [!j . An interesting reformu- 
lation of this work was recently given by Detweiler and 
Whiting [ljj, wn0 showed that the perturbation can be 
uniquely decomposed into a symmetric-singular field h%g, 
and a regular-radiative field h~g\ the full (retarded) per- 
turbation is then h a g = h^g + h^ g . Detweiler and Whit- 
ing were able to establish that while h^g reproduces the 
singularity structure of the metric perturbation, it exerts 
no force on the point particle; the gravitational self-force 
is then produced entirely by h^ g , which is a homoge- 
neous, regular, smooth field in a neighbourhood of the 
world line. 

The MiSaTaQuWa equations of motion |tj , in the 
Detweiler- Whiting formulation take the following 

form. Let z^(t) be parametric relations that describe 
the particle's world line 7, with r denoting proper time in 
the background spacetime of the central black hole. Let 
= dz^/dr be the particle's four- velocity, normalized 
with respect to the unperturbed metric: g flI/ u IJ "u v = — 1. 
Let D /dr denote covariant differentiation along the world 
line, defined with respect to a connection compatible with 
g a g. Then the particle's equations of motion are 
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where a semicolon indicates covariant differentiation with 
respect to the background connection. The right-hand 
side of Eq. (|l.l(l is the gravitational self-acceleration of 
the point particle; multiplying by m would give the grav- 
itational self-force. Equation (|l.ll) is equivalent to the 
statement that the particle moves on a geodesic in a 
spacetime with metric g a g + /i^, but the description of 
the world line refers to the background spacetime. The 
right-hand side of Eq. (|1.1|) is of order m, and the grav- 
itational self-acceleration is therefore O(m); it vanishes 
in the test-mass limit and the motion becomes geodesic 
(in the background spacetime). 

The decomposition of h a g into singular "S" and ra- 
diative "R" fields relies on a specific choice of gauge for 
the metric perturbation, which must satisfy the Lorenz 
gauge condition 

V^^-^W) =0. (1.2) 

This choice ensures that h a p satisfies a (hyperbolic) wave 
equation and that the correct, retarded solution can be 
identified. The singularity structure of the perturbation 
near the world line can then be determined by a local 
analysis (see, for example, Ref. 0), and h^g is con- 
structed without ambiguity so that it exerts no force 
on the particle. The regular field h^ 3 is then the dif- 
ference between the retarded solution and the locally- 
constructed singular field; this satisfies a homogeneous 
version of the wave equation satisfied by the full metric 
perturbation, and the metric g a p + h^g is a solution to 
the linearized Einstein field equations in vacuum. 

The Lorenz gauge therefore presents itself as a pre- 
ferred gauge for this problem, and it has been shown 
that in the first post-Newtonian approximation, Eq. 
agrees with the standard Einstein-Infeld-Hoffmann equa- 
tions of motion in a common domain of validity [1 If . But 
it is important to note that the equations of motion of 
Eq. are not gauge invariant 0: different gauge 

choices will lead to different results. 



Self-acceleration by mode sums 

A concrete evaluation of Eq. Hl.lfl is challenging and 
involves a large number of steps; for this discussion we 
consider the case of a particle orbiting a Kerr black hole. 

The first sequence of steps are concerned with the com- 
putation of the metric perturbation h a g produced by 
a point particle moving on a specified geodesic of the 
Kerr spacetime. A method for doing this was elaborated 
by Lousto and Whiting and Ori 0], buildin g on 
the pionee ring work of Teukolsky ^5|, Chrzanowski |16| . 
and Wald [l7j. The procedure consists of (i) solving the 
Teukolsky equation for one of the Newman-Penrose quan- 
tities V'o and ^4 (complex components of the Weyl tensor) 
produced by the point particle; (ii) obtaining from ?po or 
■04 a related (Hertz) potential \& by integrating an ordi- 
nary differential equation; (iii) applying to f a number 



of differential operators to obtain the metric perturba- 
tion in a radiation gauge that differs from the Lorenz 
gauge; and (iv) performing a gauge transformation from 
the radiation gauge to the Lorenz gauge. 

It is well known that the Teukolsky equation can 
be separated when ipo or ip^ is expressed as a multi- 
pole expansion, summing over modes with (spheroidal- 
harmonic) indices I and m. In fact, the procedure out- 
lined above relies heavily on this mode decomposition, 
and the metric perturbation returned at the end of the 
procedure is also expressed as a sum over modes h l a g. 
(For each I, m ranges from —I to I, and summation of m 
over this range is henceforth understood.) From these, 
mode contributions to the self-acceleration can be com- 
puted: a M [/i/] is obtained from Eq. (|1.1|) by substituting 
h l a g in place of h^g. These mode contributions do not 
diverge on the world line, but a^[hi\ is discontinuous at 
the radial position of the orbit. The sum over modes, on 
the other hand, does not converge, because the "bare" 
acceleration (constructed from the retarded field h a g) is 
formally infinite. 

The next sequence of steps is concerned with the reg- 
ularization of each a^lh i ] b y removing the contribution 
from h s af3 [3 E HJ ||1 Hi E| . The singular field can 
be constructed locally in a neighbourhood of the parti- 
cle, and then decomposed into modes of multipole order 
I. This gives rise to modes a M [/if] for the singular part 
of the self-acceleration; these are also finite and discon- 
tinuous, and their sum over I also diverges. But the true 
modes a^/ij ] = a M [ft-/] — a M [ft-f] of the self-acceleration 
are continuous at the radial position of the orbit, and 
their sum does converge to the particle's acceleration. (It 
should be noted that obtaining a mode decomposition of 
the singular field involves providing an extension of h^g 
on a sphere of constant radial coordinate, and then inte- 
grating over the angular coordinates. The arbitrariness 
of the extension introduces ambiguities in each a M [/if], 
but the ambiguity disappears after summing over I.) 

The gravitational self-acceleration is thus obtained by 
first computing a M [ft,/] from the metric perturbation de- 
rived from ipo or ^4, then computing the counterterms 
a M [ft-f] by mode-decomposing the singular field, and fi- 
nally summing over all a^/if 1 ] = a^[hi\ — a^lhf]. This 
procedure is lengthy and involved, and thus far it has not 
been brought to completion, except for the special case of 
a particle falling radially toward a nonrotating black hole 
|24j . In this regard it should be noted that replacing the 
central Kerr black hole by a Schwarzschild black hole sim- 
plifies the task considerably. In particular, because there 
exists a practical and well-developed formalism to de- 
scribe the metric perturbations of a Schwarzschild space- 
time [2^, |2(| I27], 13 |2i| , there is no necessity to rely on 
the Teukolsky formalism and the complicated reconstruc- 
tion of the metric variables. 
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Low multipoles — this work 

The procedure described above is not complete. The 
reason is that the metric perturbations 1v a p that can be 
recovered from ipo or ^4 do not by themselves sum up to 
the complete gravitational perturbation produced by the 
moving particle. Missing are the perturbations derived 
from the other Newman-Penrose quantities: ip\, ip2, and 
-03- While ipi and ip3 can always be set to zero by an 
appropriate choice of null tetrad, ip2 contains such im- 
portant physical information as the shifts in mass and 
angular-momentum parameters produced by the parti- 
cle [3(j . Because the mode decompositions of ipo an d 
ip4 start at I = 2, we might colloquially say that what 
is missing from the above procedure are the li l = and 
/ = 1" components of the metric perturbations. It is not 
currently known how the procedure can be completed so 
as to incorporate all components of the metric perturba- 
tions. 

In this paper we consider the contribution of these low 
multipoles (1 = and I = 1) to the gravitational self- 
acceleration. To make progress we shall take the central 
black hole to be nonrotating, and the metric of the back- 
ground spacetime to be a Schwarzschild solution. This 
simplification allows us to use the robust formalism of 
gravitational perturbations of the Schwarzschild space- 
time [H HE EH3 El, and more importantly, to de- 
fine precisely what is meant by the u l = and 1=1" 
modes of the perturbation field. In this context the as- 
sociations between the I = mode and a shift of mass 
parameter, the odd-parity I = 1 mode and a shift of 
angular-momentum parameter, and the reduction of the 
even-parity I = 1 mode to a gauge transformation, were 
first established by Zerilli [26|. These associations are 
central to our discussion, and we believe that the results 
derived here will have a direct counterpart in the case of 
a Kerr black hole: The missing metric perturbations of 
the Kerr spacetime will be equivalent to our I = and 
I = 1 perturbation modes in the limit where the black- 
hole angular momentum goes to zero. 

To keep our discussion concrete and the mathematical 
complexities to a minimum, we calculate the I = and 
/ = 1 perturbation modes for the specific case of a par- 
ticle moving on a circular orbit of radius R and angular 
velocity fl — ^jMjR?. While finding solutions to the rel- 
evant perturbation equations can be a simple task when 
adopting a simple choice of gauge (as we shall see), we 
insist here, for reasons that were listed before, that the 
I = and / = 1 perturbation modes should be calculated 
in the Lorenz gauge. This complicates the structure of 
the perturbation equations, and finding solutions is more 
challenging. We nevertheless are able to find exact ana- 
lytical solutions for the cases I = and / = 1 (odd parity) . 
For even-parity 1=1, however, we have to rely on nu- 
merical methods for exact results, and a post-Newtonian 
approximation for analytical results. 

In the remaining sections of the paper we calculate the 
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FIG. 1: Internal values of a[l](R), rescaled by a common factor 
of 2>m/R?. For R 2> M we have the following asymptotic 
behaviors: a<[Z = 0] ~ 3(m/R 2 )(M/R), a< [l = l;odd] ~ 
-4{m/R 2 ){M/R), and a < [l — l;even] ~ 3(m/R 2 ). An exact 
expression for a < [l — 0] appears in Eq. (|3.15^ below. An 
exact expression for a < [l = l;odd] appears in Eq. 14.21 . The 
values for a < [l = l;even] are obtained from Eq. I5.55|l and 
the results listed in Table I. 

contributions 

a^lhl—o], a M [/l;=l;odd parity], [/ii=l;cvcn parity] 

to the "bare" self-acceleration of a particle moving on a 
circular orbit around a Schwarzschild black hole. Our 
expressions are finite but discontinuous at the radial po- 
sition of the particle: the answers obtained when ap- 
proaching r = R from the interior (r < i?), and those 
obtained from exterior (r > R), do not match. We find 
in all cases that the contribution to the bare acceleration 
is purely radial: a^[hi] = a r [hi]5^ for all three modes con- 
sidered here. Moreover, in all cases the self-acceleration 
is conservative and does not contribute to the radiation 
reaction. 

To keep the notation simple we shall set 

a[l = 0] = a r [h l=0 ], 
a[i = l;odd] = a r [h l=1 . odd parity], (1.3) 
a[Z = l;even] = a r [h l=1 . ovcn parity]- 

We display our results for a[l](R) in two figures. In Fig. 1 
we show the results as calculated from the orbit's interior 
(r — > R~) and in Fig. 2 we show our results as calculated 
from the orbit's exterior (r — > R + ). These results do not 
have an immediate physical meaning. To produce mean- 
ing they must be included with higher-multipole contri- 
butions in a sum over all modes. Because there exists no 
procedure to uniquely remove the "S part" of the I = 
and I = 1 perturbation modes (as was mentioned previ- 
ously), we are not able here to produce expressions for 
the low-multipole contributions to the regularized self- 
acceleration. 

Using purely analytical methods, Nakano, Sago and 
Sasaki [3l| calculated the self-acceleration to first post- 
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FIG. 2: External values of a[/](_R), rescaled by a common 
factor of m/R 2 . For R > M we have the following asymp- 
totic behaviors: a>[? = 0] ~ —m/R 2 , a < [l = l;odd] ~ 
2{m/R 2 )(M/R), and a<[Z = l;even] ~ -3P(m/R 2 ){M/R), 
where /3 ~ 2 is numerically estimated at the end of Sec. V. 
An exact expression for a> [I — 0] appears in Eq. lEHgl below. 
An exact expression for a>[7 = l;odd] appears in Eq. I|4.3^ . 
The values for a>[Z = l;even] are obtained from Eq. 15.55|l 
and the results listed in Table I. 



rameter that occurs at r = R. We then calculate a[l = 0], 
the corresponding contribution to the self-acceleration. 

In Sec. IV we do the same for the 1 = 1, odd-parity 
perturbations. These are associated with the change of 
angular-momentum parameter that occurs across the or- 
bit, and they give rise to the contribution a[l — l;odd] 
to the self-acceleration. 

In Sec. V we consider the 1 = 1, even-parity gravita- 
tional perturbations, which are associated with the mo- 
tion of the central black hole around the system's center 
of mass. This calculation is considerably more involved 
than the others, because here the source of the pertur- 
bations is time dependent. Solving the vectorial wave 
equation that converts the perturbations from the Zerilli 
gauge to the Lorenz gauge requires numerical techniques, 
except when R 3> M and we can rely on approximate 
analytical methods. In this section we obtain exact nu- 
merical results for a[l — l;even], as well as approximate 
analytical results for R 3> M. 

In Sec. VI we discuss our results and offer a number of 
concluding remarks. 

II. NEWTONIAN SELF-ACCELERATION 



Newtonian order for circular orbits of the Schwarzschild 
geometry. For the even and odd parity I = 1 modes, their 
results for the contribution to the "bare" self-acceleration 
agree with ours at the 1PN level, as expected. It ap- 
pears that an extension of their methods to higher post- 
Newtonian orders might be substantially complicated by 
the difficulty caused by the even-parity I = 1 perturba- 
tions, the case for which we have to rely on numerical 
methods. For the I = mode our results for the "bare" 
self-accelerations disagree with theirs at 1PN. We believe 
that the discrepancy is caused by the implementation of 
boundary conditions at the event horizon, and we discuss 
this matter in some detail in the Appendix. 



In this section we consider a Newtonian system involv- 
ing a large mass M and a much smaller mass to. The 
position of the small mass relative to the center of mass 
is described by the vector R(t), while the position of the 
larger mass is described by pit). Taking the center of 
mass to be at the origin of the coordinate system, we 
have 

mR + Mp = 0. (2.1) 

We denote the position vector of an arbitrary field point 
by x, and r = \x\ is its distance from the center of mass. 
We shall also use R=\R\ and p = \p\. 



Organization of this paper 

In Sec. II we set the stage with a discussion of the 
gravitational self-force in Newtonian theory. This simple 
analogue to the relativistic problem sheds considerable 
light on the meaning of the self-acceleration and its de- 
composition into singular "S" and regular "R" fields. We 
show in particular that in Newtonian theory, the I = 1 
contribution to the self-acceleration is responsible for an 
important finite-mass correction to the particle's angu- 
lar velocity. We take this as a clear suggestion that in 
the relativistic problem, the low multipole contributions 
to the gravitational self-acceleration produce important 
physical effects. 

In Sec. Ill we compute, in the Lorenz gauge, the I = 
gravitational perturbations produced by a particle in a 
circular orbit around a Schwarzschild black hole. These 
perturbations are associated with the change of mass pa- 



Test-mass description 

We begin with a test-mass description of the situation, 
according to which the smaller mass moves in the grav- 
itational field of the larger mass which is placed at the 
origin of the coordinate system. The background New- 
tonian potential is 

$oO = 2.2 

r 

and the background gravitational field is 

M 

9o = -V$o = 5- (2.3) 

In this description, the smaller mass m moves according 
to d 2 R/dt 2 = go(x = R). If the motion is circular, then 
to possesses a uniform angular velocity given by 

«o 2 = ^, (2.4) 
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where R is the orbital radius. These results are in 
close analogy with a relativistic description in which the 
smaller mass is taken to move on a geodesic of the back- 
ground spacctime, in a test-mass approximation. 



Beyond the test-mass description: singular "S" and 
regular "R" perturbations of the Newtonian 
potential 

We next improve our description by incorporating the 
gravitational effects produced by the smaller mass. The 
exact Newtonian potential is 



M 



RV 



(2.5) 



and for m -C M this can be expressed as $(a?) = $>o(x) - 
6<$>(x), with a perturbation given by 



5$(x) 



M 



\x - p\ 



r 



\x-R\ 



(2.6) 



This gives rise to a field perturbation Sg that exerts a 
force on the smaller mass. This is the particle's "bare" 
self-acceleration, and the correspondence with the rela- 
tivistic problem is clear. 

An examination of Eq. (|2.6|) reveals that the last term 
on the right-hand side diverges at the position of the 
smaller mass. But since the gravitational field produced 
by this term is isotropic around R(t), we know that this 
field will exert no force on the particle. We conclude that 
the last term can be identified with the singular "S" part 
of the perturbation, 



$ s (x) 



las - R\ 



(2.7) 



and that the remainder makes up the regular "R" field, 

M M 
\x — p\ r 



(2.8) 



The full perturbation is then given by S&(x) = $s( a; ) + 
$n(x), and only the "R potential" affects the motion of 
the smaller mass. Once more the correspondence with 
the relativistic problem is clear. 

It is easy to check that to first order in m/M, Eq. (|2.8|) 
simplifies to 



$n(x) = m 



R x 



(2.9) 



this simplification occurs because thanks to Eq. Q2.ip . p 
is formally of order m/M <C 1. The regular "R" part of 
the field perturbation is then 



9r(x) = m 



3{R-x)x-r 2 R 



(2.10) 



and evaluating this at the particle's position yields a cor- 
rection to the background field go(x = R) = —MR/R 3 



given by gn{x = R) = 2mR/R 3 ; the force still points in 
the radial direction but the active mass has been shifted 
from M to M — 2m. For circular motion the angular 
velocity becomes 



n 2 = 



M - 2m 



R 3 



(2.11) 



This can be cast in a more recognizable form if we express 
the angular velocity in terms of the total separation s = 
R + p = (1 + m/M)R between the two masses. To first 
order in m/M we obtain fi 2 = (M + m)/s 3 , which is just 
the usual form of Kepler's third law. The regular part 
of the field perturbation is therefore responsible for the 
finite-mass correction to the angular velocity. 



Multipole decomposition of the perturbations 

We now examine the low-multipole content of SQ(x), 
$s(as), and $ R (as). 

It is evident from Eq. (|2.9|l that "^.(cc) possesses a pure 
dipolar form, and its multipole decomposition therefore 
involves a single term at I = 1. As we have seen, this 
dipole potential is responsible for an important finite- 
mass correction to the orbital frequency. We take this 
as a clear indication that in the relativistic context, the 
I = 1 contribution to the metric perturbations produces 
important physical effects that should not be ignored. 
Since there is no analogue in Newtonian theory to the 
odd parity metric perturbations, this statement might 
be restricted to the 1 = 1, even parity perturbations of 
the Schwarzschild spacetime. 

While ( E ) R(a;) possesses only a dipole component, the 
same is not true of < l ) s( a; ) and 6&(x). Their monopole 
components are given by 



1=0, 



x) 



-m/R 

—m/r 



r <R, 



r>R, 

and this gives rise to a monopole field perturbation 



(2.12) 



Sg 



1=0 







-mx/r 3 



r < R, 
r > R. 



(2.13) 



This is discontinuous at x — R: the field is zero when 
the limit is taken from the inside, and equal to —mR/R 3 
when taken from the outside. The jump in the monopole 
field perturbation is given by 



[5g l =°] = 6g l =°(x = R) 



outside 



-Sg" 



R) 



m 

R^ R - 



(2.14) 



These results, which could be described as the Newtonian 
"bare" self-acceleration for I = 0, will be recovered as 
limits of our exact relativistic expressions in Sec. III. 

The dipole component of the singular potential is cal- 
culated to be 



-m(R-x)/R 3 
-m(R ■ x)/r 3 



r<R, 
r>R, 



(2.15) 
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and adding this to Eq. (|2.9() we find 

1 1 



5<f> l=i (x) = m(R-x) 



i? 3 



(r < R) (2.16) 



and S^ l ~ 1 (x) = for r > R. This gives rise to the field 
perturbation 



Sg 1-1 = m 



3(R-x)x-r 2 R R 
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(r<R) (2.17) 



and Sg l=1 = for r > R. The dipole field also is discon- 
tinuous at x = R. it is zero when the limit is taken from 
the outside, and equal to 3mR/R 3 when taken from the 
inside. Its jump, defined as in Eq. (|2.14|> . is given by 



3m 



(2.18) 



These results, which could be described as the Newtonian 
"bare" self-acceleration for I — I, will also be recovered 
as limits of our exact relativistic expressions in Sec. V. 



III. MONOPOLE GRAVITATIONAL 
PERTURBATIONS 

Our task in this section is to calculate the I = metric 
perturbations of the Schwarzschild spacetime produced 
by a particle of mass m in circular orbit at a radius R. Wc 
shall also calculate the associated contribution to the self- 
acceleration, a[l = 0] as defined by Eq. I|1.3|) . We need 
the perturbations in the Lorenz gauge, and our strategy 
will be to obtain them first in the simpler Zerilli gauge, 
and then look for a transformation to the Lorenz gauge. 



Perturbations in the Zerilli gauge 

The monopole perturbations produced by a point par- 
ticle in arbitrary motion around a Schwarzschild black 
hole were first computed by Zerilli [26j . With his specific 
choice of gauge for circular motion, the metric perturba- 
tions are 

ft-^-jr^eir-R) (3.1) 



and 



h z rr = ^§@{r-R), (3.2) 

where / = l-2M/r, E = (1-2M/R)(1-3M/R)- 1 / 2 is 
the particle's energy per unit rest mass, and Q(r — R) is 
the Heaviside step function. It is easy to check that for 
r > R, g a p + Jt^p is another Schwarzschild metric with 
mass parameter M + mE. The perturbation therefore 
describes the sudden shift in mass parameter that occurs 
at r = R. 



Transformation to the Lorenz gauge 

The metric perturbation of Eqs. i|3.1[l and (|3.2[1 does 
not satisfy the Lorenz gauge condition of Eq. I|1.2|l . Wc 
therefore seek a vector field £ a that generates a transfor- 
mation from the Zerilli gauge to the Lorenz gauge. This 
vector must possess only an I = component, and so it 
must be of the form £ Q = [0, £(r), 0, 0]. As the pertur- 
bation is static, there is no need to include a component 
in the time direction (this point is elaborated in the Ap- 
pendix), nor a time dependence in the radial component. 

To find this vector we express the Lorenz-gauge metric 
perturbations in the standard Regge- Wheeler [2J| form 



h-AB 



= fH (r), 

= H 2 (r)/f, 
= r 2 n AB K{r), 



(3.3) 



where the upper-case latin indices run over the angular 
coordinates 9 and (j), and £Iab — diag(l, sin 2 9) is the 
metric of the unit two-sphere. We have set h tr = Hi(r) = 
on the grounds that the perturbation must be static. 
A similar notation can be used to express the Zerilli- 



gauge perturbations, and we have K z 
;ero. The gauge transfo 
2£(a;/3) and this translates to 



while Hn and 



i?2 are nonzero. The gauge transformation is given by 



H - H% + 

2/ 

r 



(3.4) 



where a prime indicates differentiation with respect to 
r. The new perturbation will satisfy the Lorenz gauge 
condition if 

,/ , , a 2M 2(2r-3M) 4/ 

(3.5) 

Using Eqs. I|3.1|) . I|3.2I) . and Q3.4[l. this becomes an ordi- 
nary differential equation for £(r): 



mE 



R-2M^- R ^ 
2mE R - 3M 



- 2 / R-2M 



Q(r - R). (3.6) 



Our task is now to find a solution to this equation. 

The function £(r) can be expressed as a superposition 
of interior and exterior solutions, 



Z(r)=U(r)e(R-r)+£ > (r)G(r-R). 



(3.7) 



The interior solution (r) satisfies the homogeneous ver- 
sion of Eq. I|3.6|l . while the exterior solution £>(r) satisfies 
Eq. i|3.6fl with S(r — R) set equal to zero and Q(r — R) 
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set equal to 1. The solutions must comply with the jump 
conditions 



iER 



(R - 2Mf 



(3.8) 



where [£] = £> (r = R) — £< (r = i?) and a similar def- 
inition holds for [£']. Equations (|3.8|1 and (|3.4|1 imply 
that in the Lorenz gauge, the metric perturbations are 
continuous at r = R: [Ho] = [H2] = [K] = 0. 

The interior solution is a linear superposition of the 
two independent solutions £1 = [r(r — 2Af)]~ 1 and £ 2 — 
r 2 /(r — 2M). Regularity at the event horizon requires 



that £ be well behaved in the limit r 
therefore choose 



U(r) 



2Mr + AM 2 



2M. We must 



(3.9) 



where a is a constant that will be determined by the jump 
conditions. The exterior solution is a linear superposition 
of £1, £2, and the particular solution £ p = —mE{R — 
3M)(R - 2M)- 1 T{r), where 

r(r) = -{[9Mr(r - 2M)] -1 [3r 3 ln(l - 2M/r) - 3Mr 2 

- \2M 2 r + 44M 3 - 24Af 3 ln[r/(2M) - l] j. (3.10) 

Because £1 ~ 1/r 2 , £2 ~ r, and r ~ 1 when r — > 00, 
proper asymptotic behavior requires that we discard £2 
from the exterior solution. We then have 

£>(r) = 6 M ,^ -^ !" 3 » r W. (3-11) 



r(r - 2M) 



R-2M 



where b is a constant that will be determined by the jump 
conditions. 

The gauge vector is now fully determined: The interior 
solution is given by Eq. (|3.9|l and the exterior solution by 
Eq. I|3.11|l with the function T(r) displayed in Eq. I3.10|l . 
The complete gauge vector field is then constructed as 
in Eq. Q3.7|l . and the constants a and b are determined 
by the jump conditions of Eq. I|3.8|) . This is sufficient 
information to calculate the Lorenz-gauge metric pertur- 
bations with the help of Eqs. I|3.3|) and (|3.4|) . Because the 
resulting expressions are moderately lengthy, we shall not 
display these results here, but proceed instead with the 
calculation of the self-acceleration. 

Before moving on we wish to call attention to the fact 
that in the foregoing manipulations, the requirements of 
staticity, regularity at the event horizon, and regularity 
at infinity have allowed us to construct a unique solution 
to the perturbation equations in the Lorenz gauge. This 
conclusion is elaborated in the Appendix. 



Monopole contribution to the self-acceleration 

The self-acceleration produced by the I = perturba- 
tions can be expressed as a sum of two terms, 

a[l = 0] = a[l = 0; Zerilli] + a[l = 0; gauge], (3.12) 



where a[l = 0; Zerilli] is the radial component of the accel- 
eration vector constructed as in Eq. 1)1. 1|) but by replacing 
w ith ^a/3> wn ile a[l = 0; gauge] is constructed from 
^gaugo _ _2^ a .^. The calculation involves the particle's 

velocity vector u>* = (1 - 3M/R)~ l / 2 [1, 0, 0, fl], and at 
the end h a p- n must be evaluated at the position of the 
particle (r = R, 9 = ir/2, and <f> = Sit), either from the 
orbit's interior (r < R) or from its exterior (r > R). This 
leads to two different values for the acceleration, a< and 
a>, respectively. Such a discontinuity was encountered 
before in a Newtonian context — refer back to Eq. H2.13f) . 
The external value of the Zerilli acceleration is given 

by 



a>[l = 0; Zerilli] 



mE 



R(R - 3M) ' 



(3.13) 



while the internal value is zero: a<[Z = 0; Zerilli] = 0. 
The gauge acceleration, on the other hand, is found to 
be 



a[l = 0; gauge] 



3M(R - 2Mf 
' R 4 {R-3M) 



and by virtue of Eq. I|3.8|l . the internal and external val- 
ues are equal. The gauge vector can most simply be eval- 
uated from the orbit's interior, and Eq. I|3.9|l gives £(i?) = 
a(R 2 + 2M R + 4M 2 ) /R. But the jump conditions imply 
a = i(mE/M)[(R-3M)ln(l-2M/R)-M]/(R-2M), 
and altogether we obtain 



a[l = Q; gauge] 



~(R-2M)(R 2 
mE 



2MR + AM 2 



M 



R - 3M 



R 5 
In I 1 



2M\ 



.(3.14) 



From Eqs. (|3.12|) - (|3.14|l we arrive at our final results. 
The internal value for the I = self-acceleration is 



a< [1 = 0] = mE 



~(R- 2M){R 2 + 2MR + AM 2 



M 



R-3M 
while the external value is 

a>[; = 0] = -mE R "- MRS 



R 5 
In I 1 



2M\ 
~R~ 



(3.15) 



8M 4 



R 5 (R-3M) 

~{R-2M)(R 2 -\ 
mE 



2MR + AM 2 ) 



R 5 



x In 1 



2M\ 



(3.16) 



When R^> M these expressions simplify to a< [I = 0] 



3mM/R 3 and a>[l = 0] m/R 2 



iM/(2R 3 ); the 



internal value is smaller than the external value by a 
factor of order M/R <C 1. These limiting expressions 
are compatible with the Newtonian results displayed in 
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Eq. (|2.13() . They differ, however, from the results of 
Nakano, Sago, and Sasaki [3l]], which are displayed in 
their Eq. (E19) — our expressions are smaller than theirs 
by a term 4mM / R 3 . This discrepancy is explained in the 
Appendix. Equations l|3.15|l and (|3.16l) were used to gen- 
erate the curves shown in Figs. 1 and 2. 



IV. DIPOLE, ODD-PARITY GRAVITATIONAL 
PERTURBATIONS 



In this section we calculate the 1 = 1, odd-parity metric 
perturbations of the Schwarzschild spacetime produced 
by a particle of mass m in circular orbit at a radius R. 
From these we shall derive their contribution to the self- 
acceleration, a[l = l;odd] as defined by Eq. (|1.3|) . Here 
we shall find that the expressions provided by Zerilli al- 
ready satisfy the Lorenz gauge condition. 

The dipole, odd-parity perturbations produced by a 
point particle in arbitrary motion around a Schwarzschild 
black hole were first computed by Zerilli [2(J and 
shown to be intimately related to the shift in angular- 
momentum parameter that occurs at the orbit. After 
specializing to circular motion in the equatorial plane, 
his results read 



httf, = —2mLsm 



r 2 /R 3 
1/r 



r < i?, 
r> R, 



(4.1) 



where L = [MR/(1 - 3M/R)] 1 / 2 is the particle's an- 
gular momentum per unit rest mass. For r < R, the 
metric g a p + h a p differs from g a p only by a gauge trans- 
formation — it is also a Schwarzschild metric with mass 
parameter M. For r > R, g a fj + h a p is a Kerr metric lin- 
earized with respect to the angular-momentum parame- 
ter a = [m/M)L. The perturbation therefore describes 
the sudden shift in angular momentum that occurs at 
r = R. 

It is easy to check that the perturbation of Eq. (|4.1() 
satisfies the Lorenz gauge condition of Eq. (|1.2f> . It is also 
easy to show that a (time-independent) gauge transfor- 
mation within the class of Lorenz gauges would produce 
a pathological behavior of the perturbation at the event 
horizon. Equation (|4.1|l therefore gives us a unique solu- 
tion to the perturbation equations in the Lorenz gauge. 

A straightforward calculation then reveals that the in- 
ternal value of the I — I, odd-parity contribution to the 
self-acceleration is 



a<[Z = 1; odd] 



4mM 1 - 2M/R 
R 3 (1-3M/R) 3 / 2 ' 

while the external value is 

2mM 1 - 2M/R 



a > [l = 1; odd] 



R 3 (1-3M/R) 3 / 2 ' 



(4.2) 



(4.3) 



These results have no analogue in Newtonian theory. 
Equations (|4.2|) and (|4.3|) were used to generate the 
curves shown in Figs. 1 and 2. 



V. DIPOLE, EVEN-PARITY GRAVITATIONAL 
PERTURBATIONS 



Our task in this section is to calculate the 1 = 1, even- 
parity metric perturbations of the Schwarzschild space- 
time produced by a particle of mass m in circular orbit at 
a radius R. We shall also calculate the associated contri- 
bution to the self-acceleration, a[l = l;even] as defined 
by Eq. (|1.3fl . Once more we need the perturbations in 
the Lorenz gauge, and as in Sec. Ill our strategy will be 
to obtain them first in the simpler Zerilli gauge, and then 
look for a transformation to the Lorenz gauge. The so- 
lution to the wave equation satisfied by the gauge vector 
field will be obtained numerically and provided in tab- 
ulated form. It will also be obtained analytically in a 
post-Newtonian expansion in powers of AI/R. 



Perturbations in the Zerilli gauge 

The dipole, even-parity perturbations produced by a 
point particle in arbitrary motion around a Schwarzschild 
black hole were first computed by Zerilli 26] in a simple 
choice of gauge. After specializing to circular motion, his 
results become 

^ R — 2M , q„o , \ 
■(l-r 3 n 2 /M) 



2mE ■ 



r(r - 2M) 
xsm6cos(4>-flt)e(r-R), (5.1) 

-6mEQ(R-2M)- 



' (r - 2M) 2 
x sin sin((/> - fit) 9(r — R), 

h% = 6mE(R-2M) ir _ r 2M)3 

x sin6»cos(</> - flt)e{r - R), 



(5.2) 



(5.3) 



where 57 = y/M/R 3 is the particle's angular velocity and 
E = (1 - 2M/R)(1 - 3M/R)- 1 / 2 is its energy per unit 
mass. Here we see that the perturbations are time depen- 
dent, and this complicates considerably the task of find- 
ing the transformation to the Lorenz gauge. Equation 
(|5.1|) reveals that the Zerilli gauge is not asymptotically 
flat, since hf t grows linearly with r as r — > oo. This indi- 
cates the fact that the metric g a p + h^g is expressed in a 
noninertial coordinate system anchored to the black hole 
instead of the system's center of mass. This statement 
will be elaborated below. 



Perturbations in a singular gauge 

The metric perturbations of Eqs. (|5.1|) - (|5.3|) do not 
satisfy the Lorenz gauge condition of Eq. (|1.2|) . To trans- 
form to the Lorenz gauge we proceed in two steps. We 
shall first transform to a gauge in which the perturbation 
is zero everywhere, except at r = R where it is singular. 
We shall then go from this singular gauge to the Lorenz 
gauge. 
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It is well known from Zerilli's work [26j that in vacuum, 
a dipole, even-parity perturbation can be completely re- 
moved by a gauge transformation. Such a perturbation, 
therefore, represents a coordinate transformation; and as 
we have already suggested, for r > R the metric g a f3+h^p 
is just a Schwarzschild solution expressed in a noninertial 
coordinate system. The perturbations of Eqs. (|5.1|I - H5.3|I . 
however, are not pure gauge because of the presence of 
the particle. They can be removed in the vacuum region 
outside of r = R, but the gauge transformation leaves 
something behind at r = R. The result of this transfor- 
mation is h s a p, the metric perturbation in what we shall 
call the singular gauge. 

The gauge transformation that removes a dipole, even- 
parity perturbation in vacuum was constructed by Zerilli 
|26j. It is generated by a vector field e a , so that h s a/3 = 



h z 



2 £ (a-,f3)- For circular motion this is given by 

et = — MR - 2M) ^ 

M v 'r-2M 

x sin0sin(0 - fit) 9(r - R), 

mE , , r 2 
e r = (R-2M)- -77 

M V '(r-2M) 2 

x sin 9 cos(0 - at) 6(r - R) , 

mE , „ r 2 

s e = (R-2M) 

6 M 'r-2M 

x cos6>cos(</» - Qt)G(r - R), 
mE . „ . r 2 



2M 

x sin0sin(0 - at) 6(r - R). 
The new metric perturbation is then 



(5.4) 
(5.5) 
(5.6) 
(5.7) 



Kr = 



flR 2 sin^sin( 
R 2 . n 



l rrl> 



mE 
~M 
^mE 
~W R - 2M 
mE 
~W 
mE 



-CH)6(r-R), (5.8) 
ys((f>- at) S(r-R), (5.9) 



R 2 cos cos(0 - at) 5(r-R), 



(5.10) 



R 2 



at)S{r-R), (5.11) 



and we see that in the singular gauge, the metric per- 
turbation is proportional to 5(r — R), which is produced 
by differentiation of the step function in e a . The gauge 
transformation therefore makes the perturbation zero 
everywhere in the vacuum region outside (and inside) 
r = R, but it contributes a singular term at the orbit. 
This illustrates the fact that the presence of matter pre- 
vents the Zerilli-gauge metric perturbation from being 
pure gauge. 



Interpretation of the gauge transformation 

The preceding discussion on coordinate systems can 
be clarified if we examine the asymptotic behavior of the 



gauge vector field in the limit r —>■ oo. In this limit 
we can seek a Newtonian interpretation of the results, 
and we shall see that in the original Zerilli gauge, the 
perturbed metric is that of a moving black hole. The 
following is patterned after a similar discussion produced 
by Zerilli |2g. 

The vector e a becomes asymptotically equal to b a in 
the limit r — > oo, where 



If 
b° 



- — (R- 2M)-^rsin(9cos(</> - Qt), 

M at 

— (R — 2M) sin 9 cos(</> - fit) , 

--T7-{R — 2M)- cos<9cos(<?(> - CH), 
M r 



mE , „ 1 
— (R-2M)^— 
M 'rwa.6 



sin(</> — Ot). 



If we introduce asymptotic Cartesian coordinates x 
rsin#cos0, y = rsin#sin</>, and z — rcos9, we have 



b l 
If 



~W^ R " 2M )^( XC0S ^ + y sin fit), 
~^-{R-2M) cos fit, 



W = - — (R-2M)sinat, 
M ' 

b z = 0. 

To give a Newtonian interpretation to these results, let 
x be the Zerilli coordinates of an arbitrary field point, 
let R(t) = (R - 2M) (cos fit, sin fit, 0) be the position 
vector of the orbiting particle, and express the preceding 
equations as 



b(t) 



"IT*®' 



&' = x ■ b(t), 



where an overdot indicates differentiation with respect to 
t. The coordinate transformation generated by b a is then 

x ncw = x + b(t), t ncw = t + x-b(t). 

We can now explain that this transformation represents a 
translation from a noninertial reference frame attached to 
the black hole to an inertial frame attached to the center 
of mass. Please refer back to Sec. II for a definition of 
the notation employed here. 

In the center-of-mass frame, the particle moves on 
a trajectory R cm (t), and the black hole moves on a 
trajectory p cm (t). In the black- hole frame we have 
Rbh(t) = R(t) and Pbh(t) = 0. The center-of-mass con- 
dition is (mE)R cm + Mp cm = 0, so we have p cm = 
— (mE/M)R cm . We also have R = R cm — p cm = 
(1 + mE/M)R cla , so that R cm = R + 0(m/M) and 
p cm = -(mE/M)R + 0(m 2 /M 2 ), or 

p cm (t) - b(t) + 0(m 2 /M 2 ). 
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The vector b(t) is therefore the position of the black hole 
relative to the center of mass, and the coordinate trans- 
formation is truly a translation from the moving frame of 
the black hole to the fixed reference frame of the center 
of mass. 



Transformation to the Lorenz gauge 

We now return to the task of transforming the metric 
perturbation from the singular gauge of Eqs. (|5.8(l - (|5.11(l 
to the Lorenz gauge. The gauge transformation is gener- 
ated by a vector field £ Q , such that 

h a p = h s aj3 - 2£( Q;/3 ) (5-12) 

is the Lorenz-gauge metric perturbation. For this to com- 
ply with Eq. 1)1.2(1 , the vector field must satisfy the inho- 
mogeneous wave equation 



□r = s a , 

where □ = is the wave operator and 

S a = Vp(hf-\g^g^h s l5 ) 
is the source term. This is given explicitly by 
ml QR 



(5.13) 



(5.14) 



S< = 



(3i? - 2M)S(r - R) 



S r 



M R-2M 

+ R(R-2M)5'(r - R) sin0sin(<?!> - fit), (5.15) 
mE 1 



M R 



S° = 



(2R - 5M)S(r - R) 

+ R(R - 2M)5'{r - R) sin 6 cos(</> - fit), (5.16) 
mE 1 

R(R - 2M)5'(r - R) cos(9cos(</> - fit), (5. 17) 



(3i? - 8M)S(r - R) 



S 4> = _ 



mE 



1 



M R 2 sm 2 
R(R - 2M)5'(r - R) 



(3i? - 8M)S(r - R) 

- nt). (5.18) 



smfsm 



To arrive at these results we have invoked the distribu- 
tional identity g(r)6'(r-R) = g(R)S'(r-R)-g'(R)S(r- 
R), where g(r) is any test function and g'(r) its derivative 
with respect to r. 

To solve Eq. I|5.13|l we decompose the vector £ Q in even- 
parity spherical harmonics of degree 1 = 1. The form of 
the source term indicates that only terms with m = ±1 
are needed, and we let 



Ut,r,9 A ) = "££(^¥±(8*), (5.19) 



U(t,r,6 A ) = ^T^r)^}^)- (5-20) 
± 

Here, the lower-case latin index a refers to the t and r 
components of the vector field, while the upper-case in- 
dex A refers to the angular components; we have set 9 A = 
(B,4>) and Y ± {8 A ) = Y™= ±l {6 A ) = T ^3/(8tt) smOe***. 

The vector S a can be decomposed in a similar way, 
and to simplify the form of the reduced wave equation 
we define the functions A ± (r), i3 ± (r), and C ± (r) by the 
relations 





1 

~2 




1 

= T 2 




1 

= T 2 



l8nmE 



iSlR^le^™ 



3 M r-2M 



8n mE 



3 M 



, (5.21) 
(5.22) 
(5.23) 



With these definitions Eq. 15.13fl becomes the following 
set of ordinary differential equations: 



J 



d 2 A ± 
dr 2 



fl 2 r 2 2 

(r - 2M) 2 ~ r(r - 2M) 



A ± - 



2M/R 



B = 



2R 2 



(r - 2M) 2 R-2M 



■S(r-R) + R 2 S'{r-R), (5.24) 



2 D± 



d 2 B 



fl 2 r 2 4(r - M) 



dr 2 [(r-2M) 2 r 2 (r-2M) 



B± + ™ 2MR „ A± + ±C± = S(r - R) + R^(r - R), (5.25) 



(r - 2M) 2 r 2 



R - 2M 



d 2 C± 2M dC ± 

dr 2 + r(r-2M) dr 



+ 



2„2 



n 2 r 



(r - 2M) 2 r(r - 2M) 



r{r - 2M) 



B = R5(r ~ R) + R 5'{r — R). (5.26) 



The steps required to obtain the Lorenz-gauge metric perturbation are therefore these: First, solve Eas. (|5.24|) - 
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(|5.26(l for the functions A ± (r), i? ± (r), and C (r); sec- 
ond, insert the solutions into Eqs. (|5.21fi - (|5.23fl . and 
these into Eqs. (f5.19fl and (f5.20f) . to construct the gauge 
vector field £ Q ; third, compute h a p using Eq. 15.12fl . 



Jump conditions and asymptotic behavior 



The solutions to Eqs. (|5.24f) - (|5.26fl can be expressed 



dB ± 
dr 

dC ± 
dr 



R(R - M) 
R- 2M ' 

R(R - AM) 
R-2M : 



(5.32) 
(5.33) 



where [ip] = ip > (r = R) — (r — R) for any function %j) 
of the radial coordinate. 

Near the event horizon the interior functions can be 
expanded as 



A ± {r) 
B ± (r) 



= A±(r)e(R - r) + A±(r)e(r - -R),(5.27) 
= B±(r)G(R - r) + £±(r)6(r - i?),(5.28) 
= C±(r)Q(R-r) + C±(r)e[r-R),(5.2Q) 



where the interior and exterior solutions satisfy the cor- 
responding homogeneous equations. To account for the 
source terms, these functions must comply with the jump 
conditions 



and 



[A±] = = [C±] = R 2 



dA ± 
dr 



2R 2 



R-2M 



(5.30) 



(5.31) 



A±(r) 



B±(r) 



C±(r) 



e Ti!ir *£a n (r-2M) n , 

oo 

n=0 
oo 

e Ti!ir *^c n (r-2M) n . 



n=0 



(5.34) 
(5.35) 
(5.36) 



These forms ensure that the vector £ Q satisfies ingoing- 
wave boundary conditions at the horizon, a necessary 
condition to obtain a retarded solution to Eq. I|5.13|l . Sub- 
stitution into Eqs. H5.24f] - I|5.26fl provides recurrence re- 
lations consisting of three coupled expressions for the a n , 
6„, and c„: 



[2Mn(n - 1) T 4if2A/ 2 (2n - l)]a„ - 4A/ 2 iT 1 6„ 

= -[n(n - 3) =F 2ifiM(4n - 5)]a„_i T 2iO(2 - n)a n _ 2 + 2MRT 1 b n - 



(5.37) 



8n 2 M 3 Ra n + [4M 2 n(n - 1) T 8ifiM 3 (2n - l)]6 n 

= -8Q 2 M 2 Ra n _ 1 - 2n 2 MRa n _ 2 - [AM{n 2 - 3n + 1) T 8iQM 2 (3n - 4)]6„_i 
- [n 2 - 5n + 2 =F 2iOM (6n - 13)]6 n _ 2 T 2«fi(3 - n)&„_ 3 - 4c„_ 2 , 



(5.38) 



2nM{n =F 4zftM)c„ = -26„_i - [n(n - 3) =F 8iQM(n - l)]c„_i =f 2i£l(n - 2)c„_ 2 . 

I 



(5.39) 



For n < 0, a„, 6„, and c n are zero. For n = and 
1, Eqs. H5.37fl - H5.39|) allow ao, Oi, and cq to be chosen 
freely. Other early coefficients in the sequences are 



bo = zkiQRao, 

7 Rao . t-j 

"1 = T U'JtQl, 



Cl 



2M 2 
Co =F iflRai 
M{1 =F 4iOM) ' 



(5.40) 
(5.41) 

(5.42) 



Similarly, for large r the exterior functions can be ex- 
panded as 



A±(r) 



ri=0 



(5.43) 



OO 

C±( r ) = e ±mr 'J2tnr- n , 



n=0 



(5.44) 
(5.45) 



These forms ensure that the vector £ a satisfies outgoing- 
wave boundary conditions at infinity, another necessary 
condition to obtain a retarded solution to Eq. I|5.13[) . Sub- 
stitution into Eqs. H5.24fl - H5.26fl provides recurrence re- 
lations consisting of three coupled expressions for the a„, 
6„, and c„: 
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±2iflna n = [(n - 2)(n + 1) ± 2ittM(2n - 3)]o n -i - 4M(n 2 - 3n + l)a„_ 2 

+ 4M 2 (n - 2)(n - 3)a„_ 3 ~ 2Mir 1 5„_i, (5.46) 

± 2ifin6„ = 2n 2 RMa n - 1 + [n 2 - n - 4 ± 2ifiM (2n - 3)]S n _ x - 4Af (n 2 - 3n - l)S n _ 2 

+ 4Af 2 (n 2 - 5n + 4)6„_ 3 + 4c„_i - 16Mc„_ 2 + 16M 2 c m _ 3 , (5.47) 

and 

±2innc n = 26„_i + (n - 2)(n + l)c„_i - 2Mn(n - 2)c„_ 2 - (5.48) 



For n < 0, a„, &„, and c„ are zero. Eqs. (|5.46(l - (|5.48|l 
allow ao, bo, and Co to be chosen freely. Other early 
coefficients in the sequences are given by 

±ifld 1 = -(l±mM)a -Mb /R, (5.49) 
±iflbi = -MQ 2 Ra - (2±ifiM)6 + 2c ,(5.50) 
±if2ci = &o ~ Co- (5.51) 

The set of homogeneous solutions to Eqs. H5.24fl - 15.26fl . 
inside and outside the orbit, forms a six dimensional lin- 
ear vector space. The six amplitudes 

o-o, a i, Co, ao, bo, Co 

determine one complete homogeneous solution and may 
be considered to be the "components" of any member of 
this vector space. The six amplitudes that generate the 
particular solution which satisfies the matching condi- 
tions of Eqs. H5.3Ufl ~ H5.33fl therefore identify the member 
of the vector space that corresponds to the desired gauge 
transformation. 

Numerical integration of the ABC equations 

The numerical integration of the homogeneous ver- 
sions of Eqs. (f5.24fl - (f5.26fl is performed by first choosing 
starting points r m ; n and r max which are close enough to 
their limiting values, respectively 2M and infinity, that 
the expansions and recursion relations H5.34fl - lf5.39fl and 
H5.43l) - H5.48fl provide appropriate initial conditions for A, 
B, and C at machine accuracy with a reasonable number 
of terms in the sums (no more than 30 in our case) . Also, 
the starting points are chosen to be sufficiently close to R 
that the resulting integration to R takes only a few sec- 
onds of machine time. Satisfying these two requirements 
simultaneously, both inside and outside the orbit, is not 
difficult in practice. 

The integration routine requires six input parameters, 
the complex amplitudes ao, ai, cq, clq, bo, and to. These 
must be chosen so that the six jump conditions (|5 . 30|> 
H5.33fl are enforced. We have six algebraic equations for 
six unknowns. We pick a set of six linearly indepen- 
dent "basis solutions" , each of which has only one of the 



ao . . . cq equal to 1, all other amplitudes being zero. After 
integrating the basis solutions to R we collect the values 
of A, B, C, dA/dr, dB/dr and dC/dr, all evaluated at 
_R, in a matrix 





-A 2< 


-A 3< 


Ai> 


A 2> 


A 3 > 


\ 


-Bi< 


—B 2 < 


-£?3< 


Bi> 


B 2> 


Bz> 


-Ci< 


— C*2< 




Ci> 


C 2> 


c 3 > 




~A' X< 


A' 


A' 


A[> 


A' 


A' 




-B[< 


~ B 2< 


~ B 3< 


B[> 


B' 2> 


B 3> 


) 


-c x< 


~c 2 < 


~^3< 


C'i> 


c 2> 


c 3 > 



where we use an obvious notation; for example, A' 1K is 
the value of dA jdr at r = R for the first of the internal 
basis solutions. The required amplitudes of the six basis 
solutions form the unknown column vector x, and the 
column vector j contains the values of the discontinuities 
obtained from the jump conditions H5.3Ufl - H5.33fl . After 
integrating the six basis solutions, we are left to solve the 
system of linear equations 

Mx=j (5.52) 

for the desired amplitudes x of our basis solutions; 
these then combine to give us the desired solution of 
Eqs. H5.24fi - H5.26fl with appropriate boundary conditions. 

In our numerical work we use double-precision arith- 
metic and have adopted two different ODE integration 
routines from Chapter 16 of Numerical Recipes |32j, the 
Runge-Kutta and the Burlish-Stoer algorithms. Each of 
these routines contains an accuracy parameter. A com- 
parison of the numerical results over a range of values 
of this parameter allows us to be certain that all digits 
quoted in Table I are significant. We tested the consis- 
tency of the integrations versus the expansions by nu- 
merically integrating over a wide range in r where the 
expansions give accurate values for A, B and C . The 
consistency of the expansion routines with the integra- 
tion routines is strong evidence that coding errors have 
been eliminated. Furthermore, we have written two in- 
dependent codes, one per author, and all results were 
obtained independently before they were compared with 
each other. The agreement was well within the numer- 
ical errors of each code. Our final results for A^(R), 
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TABLE I: Computed values for the internal functions A~^(R), B<(R) and C^(R). The external values are obtained by applying 
the jump conditions: ReL4>] = Re[A<] + R 2 and the imaginary parts are identical (similar statements hold for B and C). The 
functions A~(R), B~(R), and C~(R) are obtained by complex conjugation. All digits provided are significant. Note that we 
have set M = 1 in our computations. 



R 


Re[A+] 


lm[A+] 


Re[B+] 


Im[B+] R 


Re[C^] 


Im[C+] 


6 


-39.427067 


-0.68518043 


-28.037347 


3.1558616 


-26.185013 


3.8814154 


7 


-52.930571 


-0.68313011 


-38.997736 


3.6225886 


-37.151392 


4.3105572 


8 


-68.389381 


-0.66104724 


-51.991922 


4.0456669 


-50.152019 


4.6962809 


9 


-85.826868 


-0.63522323 


-67.000673 


4.4373057 


-65.166204 


5.0547070 


10 


-105.25284 


-0.61006840 


-84.015966 


4.8042522 


-82.185790 


5.3926745 


11 


-126.67197 


-0.58677222 


-103.03408 


5.1508692 


-101.20722 


5.7139596 


12 


-150.08673 


-0.56553614 


-124.05317 


5.4802524 


-122.22887 


6.0210382 


13 


-175.49853 


-0.54624909 


-147.07229 


5.7947312 


-145.24995 


6.3157134 


11 


-202.90822 


-0.52871248 


-172.09095 


6.0961288 


-170.27013 


6.5993856 


15 


-232.31636 


-0.51271828 


-199.10891 


6.3859127 


-197.28925 


6.8731864 


16 


-263.72332 


-0.49807505 


-228.12606 


6.6652898 


-226.30729 


7.1380554 


17 


-297.12936 


-0.48461496 


-259.14236 


6.9352693 


-257.32427 


7.3947870 


18 


-332.53466 


-0.47219375 


-292.15781 


7.1967063 


-290.34024 


7.6440626 


19 


-369.93935 


-0.46068823 


-327.17244 


7.4503337 


-325.35525 


7.8864727 


20 


-409.34355 


-0.44999342 


-364.18629 


7.6967858 


-362.36939 


8.1225343 


25 


-636.35940 


-0.40594758 


-579.24539 


8.8386535 


-577.42896 


9.2219041 


30 


-913.37000 


-0.37286764 


-844.29131 


9.8621836 


-842.47459 


10.213754 


35 


-1240.3777 


-0.34683600 


-1159.3279 


10.797069 


-1157.5106 


11.123814 


40 


-1617.3835 


-0.32564845 


-1524.3578 


11.662460 


-1522.5399 


11.969054 


45 


-2044.3881 


-0.30796022 


-1939.3826 


12.471537 


-1937.5641 


12.761341 


50 


-2521.3918 


-0.29289834 


-2404.4036 


13.233822 


-2402.5846 


13.509349 


55 


-3048.3949 


-0.27986801 


-2919.4216 


13.956447 


-2917.6021 


14.219636 


60 


-3625.3975 


-0.26844793 


-3484.4373 


14.644917 


-3482.6173 


14.897304 


65 


-4252.3997 


-0.25833015 


-4099.4510 


15.303582 


-4097.6306 


15.546405 


70 


-4929.4016 


-0.24928356 


-4764.4632 


15.935949 


-4762.6423 


16.170225 


75 


-5656.4033 


-0.24113082 


-5479.4740 


16.544892 


-5477.6528 


16.771470 


80 


-6433.4048 


-0.23373328 


-6244.4838 


17.132801 


-6242.6622 


17.352398 


85 


-7260.4061 


-0.22698069 


-7059.4926 


17.701689 


-7057.6707 


17.914916 


90 


-8137.4073 


-0.22078417 


-7924.5006 


18.253269 


-7922.6784 


18.460652 


95 


-9064.4083 


-0.21507114 


-8839.5079 


18.789014 


-8837.6855 


18.991011 


100 


-10041.409 


-0.20978162 


-9804.5146 


19.310200 


-9802.6919 


19.507212 



B<(R) and C<(i?) are listed in Table I for selected val- 
ues of R. Results for A<(R), B<(R) and C<(i?) are 
obtained by complex conjugation. Results for A>(R), 



B±{R) and C±{R) are obtained from the jump condi- 
tions of Eq. (JEinj). 



Calculation of the self-acceleration 

Sub stitution of Eqs. (|537|) - (jF^|) into Eqs. ifOTJl 
these into Eqs. (|3TT~^f) . (|^U|) . and finally, these 
into Eq. 1)5. 12[) yields 

Kp = h a ap - {\^ a ]r tP +r ta [^])§{r-R) 



in an obvious notation; for example Q is the internal 
(r < R) solution to Eq. (|5.13|) . constructed from A < (r), 
B < (r), and C<(r). The first three terms on the right- 
hand side appear to be singular, but it is easy to check 
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that by virtue of Eqs. (|5.8|l - (|5.11() and l|5.30() . the factors 
multiplying S(r — R) are all zero. We therefore have 

= - + $. a ) Q(R-r)- (t> + 0(r - i?). 

(5.53) 

The jump conditions (|5.31() - (|5.33() also enforce 

and we see that in the Lorenz gauge, the metric perturba- 
tion is continuous at r — R. Equation (|5.53|) also reveals 
that the internal (r < R) and external (r > R) forms of 
h a /3 are obtained by a pure gauge transformation. The 
internal and external transformations, however, are dis- 
tinct, and the perturbation is not globally pure gauge. 
Differentiation of h a p gives 



h a p n = -[€a-,p + €p;a]r,*/8(r - R) 



'A 



R). 



Once more the singular terms vanish and we end up with 
the nonsingular (but discontinuous) tensor 



hap;-) 



■8 



>,« 7 / e ( r 



r) 

-R). 



(5.54) 



This can now be substituted into Eq. I|l.l|l to obtain the 
1 = 1, even-parity contribution to the self-acceleration; 
the calculation also involves the particle's velocity vec- 
tor, = (1 - 3M/iJ) _1 / 2 [l,0,0,n]. In the notation of 
Eq. |JOJ, we have 

a[l = l;even] = R *~™ Re[B+(R)], 

(5.55) 

which must be evaluated on either side of r = R. To 
arrive at Eq. I|5.55|l we have used the property that 
B~{R) is the complex conjugate of B + (R), so that 
B+(R) + B~(R) = 2Re[B+(R)]. That the acceleration 
vector depends only on the radial component of £ a is 
a consequence of the facts that the acceleration is pure 
gauge (in the sense given above) and that the motion is 
circular. The curves displayed in Figs. 1 and 2 were ob- 
tained by substituting the numerical results of Table I 
into Eq. Q5.55p . 

The I = 1, even-parity contribution to the self- 
acceleration takes different values depending on whether 
B + (R) is evaluated from inside or outside the orbit. By 
virtue of Eq. 1)5. 3U[) . its jump across the orbit is given by 



\a[l = l;even]] = —3mE 



R - 2M 
R 2 {R-3M)' 



(5.56) 



When R 3> M, this agrees with the Newtonian result of 
Eq. 



Self-acceleration in the post-Newtonian limit 

While we have not been able to find exact analytic so- 
lutions to Eqs. (|5.24|l - (|5.2()l) . it is possible to make some 
progress by linearizing the equations with respect to M. 
Solutions to these equations are then post-Newtonian 
approximations to the exact, numerically obtained so- 
lutions. We now set out to obtain these approximations, 
and to compare them with the numerical results. 

After linearization — recall that f2 2 = M/R 3 is linear 
in M — the homogeneous equations become 

£ + » i^-U^B.o. (5.57) 



d 2 B 
dr 2 



n 2 



4M\ 



B+^rC = 0, 



and 



d 2 C 2M dC 



dr 2 



dr 



9. 2 



AM 



C- 



AM 



(5.58) 



B = 0, 



(5.59) 

where we have omitted the ± labels for ease of notation. 
The jump conditions reduce to [A] = [B] = [C] = R 2 and 



dA 



dr 
dB 



dr 
dC 



( 2M 

= 2 n i+ ir 



M 

^ R[1+ R 



(5.60) 



dr 



= R 11- 



2M 



We notice that the equations for B and C decouple 
from the equation for A. In the sequel we will con- 
struct solutions to the B and C equations, and leave A(r) 
undetermined; for the purposes of calculating the self- 
acceleration, only B(r) is required. Our solutions will 
satisfy out going- wave boundary conditions at r — * oo, so 
that in the following, B(r) = B+(r) and C(r) = C+(r). 

To decouple the B and C equations we introduce the 
new dependent variables ip- = [B — (1 — M /r)C]/ 'R 2 and 
i'+ = [|B + §(1 - M/r)C]/R 2 , such that 



B = R [ ip + + -V- 



and 



(5.61) 



(5.62) 



Away from r = R, these functions satisfy the differential 
equations 



6 6^ 3 

z 2 z 



(5.63) 
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and 

+ = 0, (5.64) 

where we have introduced the rescaled independent vari- 
able z — fir and the small quantity v 3 — MQ — 
(M/i?) 3 ^ 2 ; a prime indicates differentiation with respect 
to z. In terms of the new variables, the jump conditions 
become 

[V>_] = v 2 , [ip'_] = 3v (5.65) 

and 

[^ + ]=l- 2 -v\ [<]=J(l- V 2 ); (5.66) 

matching is carried out at z = v. 

To find solutions to Eq. 1)5. 63[) we use the fact that v 
is small and write 

ip- = i/j + v 3 iJi + O(v 6 ). (5.67) 

Substitution into Eq. 1)5. 63[) yields an equation for tf} , 

^o + (l-^)lfc> = 0, (5.68) 
and another equation for ipi, 

(l-^)^i = ^o. (5.69) 

We first solve these equations in the domain z < v. 
Among all possible solutions to Eq. H5.63|) , we choose one 
which does not diverge in the limit z — > 0. While this con- 
dition seems appropriate for our purposes, it is important 
to understand that we cannot fully justify it here: this 
choice must be introduced as an additional assumption. 
The reason is as follows: Linearization of the equations 
with respect to M implies that Eqs. I|5.63[l and (|5.64|) 
apply only in the domain r > M, or z > w 3 , and this 
restriction prevents us from imposing a proper ingoing- 
wave condition at the horizon (r = 2M, or z = 2v 3 ). We 
must therefore identify a suitable replacement for this 
boundary condition, in the form of an asymptotic condi- 
tion holding when z is restricted by v 3 <C z <C v. Previ- 
ous experience [33|, |34( with solving the Regge- Wheeler 
equation (2{| in the low- frequency limit (MQ = d 3 « 1) 
suggests that an appropriate substitution is a regularity 
condition in the formal limit z — > 0. This is the choice 
we make here, without confirmation that this conclusion 
applies to the system H5.24H - H5.26fl . 

A regular solution to Eq. I|5.68|l is ipo(z) = zj 2 {z), or 

/ 3 \ 3 
(z) = f — — 1 J sin z cos z, (5.70) 

where j%{z) is a spherical Bessel function. Substituting 
Eq. H5.70fl into Eq. (|5.69|) and integrating returns a linear 
superposition of zj2(z), zn,2(z), and tp p (z), a particular 



solution to the differential equation. The term involv- 
ing zj 2 (z) can be discarded, as it simply renormalizcs 
the zeroth-order solution of Eq. (|5.70|) . The coefficient 
in front of zn 2 (z) must then be chosen so as to yield a 
regular solution. This gives 

^ (*) = -\ ( \ ~ Jr) sin z + 5(1 - J) cos ( 5 - 71 ) 

and the complete interior solution to Eq. I|5.63|l is 

iP<(z) = a[^ +v 3 ^+0{v 6 )]. (5.72) 

The amplitude a will be determined by matching. 

We next turn to the domain z > v and construct an 
exterior solution to Eq. H5.63[> : this will be required to 
satisfy an outgoing- wave condition as z — » 00. The pro- 
cedure is largely the same as for the interior solution, but 
is simplified by the fact that the outer boundary is part 
of the domain z S> v 3 . An outgoing- wave solution to 
Eq. H5.68fl is i>a{z) — —izh 2 {z), or 

^(z)=(l + j ~^)e iz - (5-73) 

Substituting this into Eq. (|5.69|) and integrating returns 

a linear superposition of zh 2 (z), zh 2 2 \z), and ip p (z), a 

particular solution to the differential equation. As before 

the term involving zh 2 ( z ) can be discarded, and the 
(2) 

term involving zh 2 (z) must also be eliminated because 
it represents an incoming wave. We are left with the 
particular solution, 

^ , „ Si ( 1 2i \ - 

^ )= T(? + ?r (5j4) 

The complete exterior solution to Eq. H5.63J) is then 

i>>{z)=b[^ +v 3 *p> + 0(v 6 )], (5.75) 

and the amplitude b will be determined by matching. 

The constants a and b are determined by insert- 
ing Eqs. (|5.72(l and l|5.75|i into the jump conditions of 
Eq. (|5.65fi . The results are moderately complicated, and 
we shall not display them here. The expressions, how- 
ever, simplify once we take into account the fact that v 
is small. At the matching point we find 

^<(v) = -v 2 + 0(v 4 ), ip>(v) = 0(v 4 ). (5.76) 

In the interior domain [z < v) we can take advantage of 
the fact that z is formally of order v to derive 

$<{g) = -{z/v) 3 v 2 + 0(v 4 ). (5.77) 

We now proceed with finding interior and exterior so- 
lutions to Eq. (|5.64|) . This is a much simpler task, but as 
we shall see, our solutions will not be fully determined. 
For an interior solution we write 

V>+ (z) = a (sin z — f3v 3 cos z) , (5.78) 
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where a and (3 are constants; the scaling of the cosine 
term with v 3 is introduced for convenience, in anticipa- 
tion of later results. We note that this solution is regular 
in the formal limit z — > 0, in agreement with the discus- 
sion given previously, and that it involves two undeter- 
mined constants. For an exterior solution we choose 



(5.79) 



where 7 is another constant. Substitution of Eqs. (|5.78|) 
and H5.79fl into the jump conditions of Eq. I|5.66[l allows 
us to determine a and 7, but (3 is left over as a free 
parameter. Once more the resulting expressions are too 
complicated to be displayed, but they simplify for »Cl. 
At the matching point we find 



(5.80) 



^+(") = -1+ 1/3 + -jv' + 0(iv J ,v' 1 ), 

^>(v) = (3v 2 + 0(iv z ,v A ), 
and in the interior domain we have 



ipl(z) = -{z/v) + 



P+\(z/v)+ l -{z/vf 



v 2 + 0{iv 3 ,v 4 ). 

(5.81) 

Su bstitu tion of Eqs. 1 (5775) 1 and ((51"0")) into Eqs. (|51HT) 
and H5.62f) yields 



£<(i?) = C*<(i?) = -R 2 



1 - (3— + 0(iMQ, M 2 /R 2 ) 



as well as B > (R) = C> (R) = (3MR + 
these results, Eq. I|5.55|l becomes 



(5.82) 
According to 



a< [I = 1; even] 



3mE 



(5.83) 



and 



These solutions are parameterized by (3, which cannot be 
determined here because of our lack of control over the 
behavior of the solutions near r = 2M. 

We can, however, estimate the value of j3 by fitting 
Eq. (|5.82l) to our numerical results. We proceed as fol- 
lows. First, we fit the expression 1 - (3M/ R + (3 B (M/ R) 2 
to our numerical values for B < (R)/{—R 2 ) in the inter- 
val 20 < R/M < 100; this yields (3 = 1.9936 ± 0.0006 
and (3' B = 4.07 ± 0.02. Second, we fit the expression 
1 - (3M/R + (3' C {M/R) 2 to our numerical values for 
C < (R)/(—R 2 ) restricted to the same interval; this yields 
(3 = 1.9936 ± 0.0006 and f3' c = 2.26 ± 0.02. Third, we 
fit the expression (3 + (3 B (M/R) to our numerical val- 
ues for B>(R)/(MR); this yields (3 = 1.9951 ± 0.0004 
and (3 B = —4.12 ± 0.02. Finally, we fit the expression 
/3 + /3'c(M/R) to our numerical values for C > (R)/(MR); 
this yields (3 = 1.9952 ± 0.0004 and = -2.31 ± 0.02. 
We notice an excellent consistency among our estimates 
of P, and we conclude that according to our numerical 
results, 

13 = 1.994 ±0.001. 

It is probable that the actual value is (3 = 2, and that the 
slight discrepancy results from a failure to include addi- 
tional terms in the expansions in powers of M/R. The 
two-parameter fits presented here were obtained with a 
nonlinear least-squares Marquardt-Levenberg algorithm, 
as implemented in the software gnuplot. 

The quality of the fits can be judged by comparing the 
numerically-obtained functions B < (r) and C<(r) with 
the post-Newtonian approximations of Eqs. (|5.85|l and 
(|5.86(l . in which we substitute (3 = 2. We present this 
comparison in Fig. 3 for R = 25M. We see that the ana- 
lytic expressions are very accurate for all values of r < R 
except near r = 2M. 



VI. DISCUSSION 



a> [I — 1 ; even] 



3mE j3M 
~R 2 R~' 



(5.84) 



This is compatible with the Newtonian results presented 

in Eq. (|2~T7)l . 

Substitution of Eqs. j577Jt and ((5in)) into Eqs. (|51HT) 
and H5.62fl gives us expressions for the interior functions: 



B<(r) 



hO(iMtt,M 2 /R 2 



(3 + \{r/R)-\{r/Rf 



M 
~R 

(5.85) 



and 



C<(r) 



R 2 \ (r/R) - 



(3-1 + -(r/R) 



^ + 0(iMn, M 2 /R 2 ) }. i r,.s(i. 



Using the tensor harmonic decomposition of Regge and 
Wheeler j2~j, Zerilli (2~| and many others have studied 
the metric perturbations resulting from the geodesic mo- 
tion of a small mass in the background geometry of a 
Schwarzschild black hole. Most of the attention was de- 
voted to the radiating modes, those with I > 2, and their 
analysis typically involves numerical work. 

Much less attention has been garnered by the non- 
radiating, / = and I = 1 modes. In fact, the vacuum 
I = 1 even-parity metric perturbations were shown to be 
just gauge by Zerilli. This mode contains no gravitational 
radiation, and is usually ignored in analyses involving 
gravitational perturbations of black holes. Nonetheless, 
this very mode plays an important role in self-force cal- 
culations: only the dipole mode has a Newtonian-order 
contribution to the self force, as we have shown in Sec. II. 
The dipole metric perturbations cannot be ignored. 

Zerilli found analytic expressions for the / = and 
I = 1 metric perturbations in a convenient gauge. The 
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FIG. 3: Accuracy of the post-Newtonian expressions for the 
functions B < (r) and C<(r), for R = 25M. The solid curve 
is a plot of Re[_B<(r)], as given by Eq. |5.85p . divided by the 
numerical results listed in Table I. The dashed curve is a plot 
of Re[C<(r)], as given by Eq. 15.861 . divided by the numerical 
results listed in Table I. In both cases we have set /3 = 2. The 
error is estimated to be of order (M/R) 2 ~ 0.002. The plots 
reveal that this estimate is accurate for all values of r except 
near r = 2M. 



Lorenz gauge, however, with its hyperbolic wave opera- 
tor, is preferred for self- force calculations. We, as well 
as Nakano, Sago, and Sasaki [3J, have found analytic 
expressions for the I = and odd-parity I = 1 cases. Our 
analysis of the even-parity I — 1 case is mostly numerical, 
but our procedure is robust and easy to implement. 

While the Lorenz-gauge treatment of these non- 
radiating modes is now in hand, this analysis is but a 
small part of a complete computation of the regularized 
self-acceleration, a program that was outlined in Sec. I. 
And the ultimate goal of incorporating the equations of 
motion, with their corrections of order m/M, into a wave- 
generation formalism to obtain accurate gravitational- 
wave templates, remains elusive. 

For example, the conservative forces discussed in this 
paper affect the trajectory of the small mass at order 
m/M. But the description of this effect inherently de- 
pends upon the choice of gauge. While the actual obser- 
vation of a gravitational-wave signal at a large distance 
from the system is a gauge-independent measurement, 
the details of the conversion from the self-force, as mea- 
sured in the Lorentz gauge, to the m/M corrections to 
the wave forms, which are gauge invariant, are not yet 
known. 
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APPENDIX A: MONOPOLE GAUGE 
TRANSFORMATION WITHIN THE LORENZ 
GAUGE 

In section III we calculated the self-acceleration from 
the monopole metric perturbation in the Lorenz gauge 
and claimed that this result was unique. Nevertheless, 
our results differ from the post-Newtonian results of 
Nakano, Sago, and Sasaki (NSS) |3lJ, who also work in 
the Lorenz gauge. In this Appendix we outline a possible 
cause for the discrepancy. 

Both groups begin with differing solutions in the Zer- 
illi gauge and then find differing gauge transformations 
to the Lorenz gauge, with resulting metric perturbations 
that yield differing accelerations. We ask: do our re- 
sults differ from NSS by a gauge transformation from one 
Lorenz gauge (ours) to another Lorenz gauge (theirs)? 
An affirmative answer would invalidate our statement 
that our choice of Lorenz gauge is unique. We shall 
instead argue that while our results are indeed related 
by a gauge transformation, this transformation takes our 
Lorenz gauge into another Lorenz gauge that fails to be 
regular on the event horizon. We are therefore correct 
in stating that our gauge choice is unique, because the 
gauge employed by NSS, while appropriate for a post- 
Newtonian treatment, does not have a proper relativistic 
generalization. 

Before we investigate this matter, we note that our 
expression for the metric perturbation in the Zerilli gauge 
has 



2mE 



1 



2M 



R - 2M 



Q(r-R), (Al) 



which is zero inside the orbit; this property simplifies 
our implementation of the boundary conditions at the 
horizon. Nakano, Sago, and Sasaki, their Eq. (E5), have 
instead 



'NSS 



2mE 



2M 



R - 2M 



Q(R 



6(r - R) , (A2) 
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which inside the orbit is very similar to the Newtonian 
result of Eq. lj2~T2j) . The difference 



h 



NSS 



hft = 



2mE r - 2M 
r R- 2M 



is a gauge transformation generated by 
with a = mE/(R - 2M), so that 
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lNSS 

'Ht 



2mE r - 2M 
r R- 2M ' 



(A3) 
[at, 0, 0, 0] 

(A4) 



All other components of £( a ;/3) are zero. 

Now we seek an answer to our earlier questions. We 
assume that the I — metric perturbation has already 
been obtained in one Lorenz gauge (ours, the results ap- 
pearing in Sec. Ill), and we ask whether it is possible to 
take it to another Lorenz gauge (which we imagine to be 
a relativistic generalization of the choice made by NSS) . 
We consider the most general I = gauge vector 



3 



[at,j(r),0,0] 



(A5) 



that keeps the metric perturbation static, where a is now 
an arbitrary constant. The gauge transformation gener- 
ated by this vector produces a shift in the metric per- 
turbation given by Ah a p — —2j^ a .py For the shifted 
perturbation to satisfy the Lorenz gauge condition, the 
gauge vector must satisfy the wave equation Dj a = 0. 
This gives 



.„ 2 ., 2 . _ 2aM 



(A6) 



where / = 1 — 2M/r. The general solution to this wave 
equation is 



c\r + c 2 - 



M 3 



a 

9^ 



3r 3 ln(l-2Af/r)-3Mr 2 



-\2M 2 r + 44M 3 - 24M 3 ln(r/2M - 1) 



(A7) 



where c\ and c 2 are arbitrary, dimensionlcss constants. 
Notice that the complicated function within the square 
brackets is closely related to the function V(r) defined 
in Eq. (|3~TU|) . namely [ ] = -9Mr 2 fT(r). The resulting 
non-zero shifts in the metric perturbation are 



Ah tt 



2M 

2a f + —j{r), 



r- 

2 dj(r) 



2M 



/ dr r 2 f 
n 2 6Ah u = -2rj(r). 



(A8) 



We must now find constants c\, C2, and a which yield a 
regular metric perturbation everywhere, including on the 
event horizon and at infinity. 

Behavior at the event horizon is easily examined in 
Eddington-Finkelstein coordinates 



V = t + r + 2Mln(r/2M - 1), 



K = r. 



(A9) 



The coordinates V and 1Z are well defined everywhere 
in the vicinity of the future horizon, which is located at 
1Z = 2M . One component of j a in Eddington-Finkelstein 
coordinates is 

J V = f+f/f 

= a[V-TZ- 2M\n(K/2M - 1)] +j(ft)(A10) 

With a substitution from Eq. (|A7|) it is seen that regu- 
larity of j v at the future horizon requires that 



1 2 
2ci + -c 2 + -a = 0, 



(All) 



and with this same condition j 11 = O(f). We still need 
to check the regularity of Ah a p in Eddington-Finkelstein 
coordinates, 

A/iyv = Ah tt , 
Ah VR = Ahull 

AhuK = Ahu/f 2 + Ah rr . (A12) 
With the choice of constants in Eq. (|A11|) , we see that 
Ah tt = af[l + \n{r/2M - 1)] + 0(f 2 ), (A13) 



and 



Ah tt /f 2 + Ah rr = 2a + 0(f). 



(A14) 



Thus A/iyv an d Ah-jm are regular on the future horizon, 
but A/iy-R. is singular if a is not zero. Further analysis 
shows that a choice of constants which does not satisfy 
Eq. IjAllfl only makes the shifts more singular. Exami- 
nation of behavior on the past horizon leads to the same 
conclusion. The only condition on the constants that 
makes Ah a p a regular tensor field on the horizon is a = 
and ci = — C2/8; but Ah a p is then ill behaved as r — > 00. 
We conclude that there is no monopole gauge transfor- 
mation that simultaneously preserves the Lorenz gauge 
condition and behaves properly on the event horizon and 
at infinity. This confirms that our claim was true: our 
choice of Lorenz gauge is indeed unique. 

Nonetheless, if we set c\ = and make a gauge 
transformation with Eq. (| A5|> . we find that the result- 
ing change in the self-acceleration is completely due, at 
first post-Newtonian order, to j = at; the contribu- 
tion from j r = j(r) appears at higher post-Newtonian 
order. With the value of a set to mE/(R — 2M), the 
first post-Newtonian contribution to the self-acceleration 
is —AmM j R? , and it completely accounts for the dif- 
ference between our results and those of Nakano, Sago, 
and Sasaki |3l| . who present first post-Newtonian re- 
sults in their Eq. (E19). From the discussion provided 
above, we conclude that our results are related by a gauge 
transformation, but that this transformation takes our 
Lorenz gauge into a Lorenz gauge that fails to be reg- 
ular on the event horizon. In other words, the Lorenz 
gauge employed by NSS, while appropriate for their post- 
Newtonian treatment, does not have a proper relativistic 
generalization. 
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